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Vibrational modes and low-temperature thermal properties 
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We present a phenomenological force-constant model developed for the description of lattice dy- 
namics of sp^ hybridized carbon networks. Within this model approach, we introduce a new set 
of parameters to calculate the phonon dispersion of graphene by fitting the ab initio dispersion. 
Vibrational modes of carbon nanotubes are obtained by folding the 2D dispersion of graphene and 
applying special corrections for the low-frequency modes. Particular attention is paid to the exact 
dispersion law of the acoustic modes, which determine the low-frequency thermal properties and 
reveal quantum size effects in carbon nanotubes. On the basis of the resulting phonon spectra, we 
calculate the specific heat and the thermal conductance for several achiral nanotubes of different 
diameter. Through the temperature dependence of the specific heat we demonstrate that phonon 
spectra of carbon nanotubes show one- dimensional behavior and that the phonon subbands are 
quantized at low temperatures. Consequently, we prove the quantization of the phonon thermal 
conductance by means of an analysis based on the Landauer theory of heat transport. 



I. INTRODUCTION 

Phonons play a fundamental role in the physics and 
the characterization of graphene and carbon nanotubes. 
Phenomena such as chargeji^i^ spini and heat trans- 
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infrared and Raman spectra 



9.10.11 



electron- 



phonon scatteringi^ii^ii^ii^ii^iiiii^ and its related effects 
as superconductivityiSi and resistivity22i2i can be under- 
stood, in most situations, only with a detailed knowledge 
of the phonon spectrum. In particular, much effort has 
been done for determining thermal and transport prop- 
erties, which closely depend on the vibrational modes. 
The most striking results of experimental and theoret- 
ical research in this domain are the observation of the 
quantization of the phonon band structure through an 
analysis of the specific heat^ the discovery of ballistic 
phonon transport^^ and the measurement of the quantum 
of thermal conductance in a nanowirei^ In technological 
applications such as nanotube-based electronic devices, 
thermal properties are of central importance for under- 
standing and controlling heat dissipation and self-heating 
effects <^ Efficient thermal management is required for 
ensuring the performance and stability of the devices. 

Much experimental work has been done for measuring 
vibrational spectra— liii^i^ and for detecting and con- 
trolling the phonon population of isolated nanotubes— 
The best known feature of experimental data is the strong 
Raman-active radial breathing mode (REM), which is 
often used for the characterization or identification of 
different nanotubes in a sample— From a theoretical 
point of view, phonon modes of graphene have been stud- 
ied either by effective models^i^^ or by ah initio calcula- 
tions— 1^'^>^>^>^'^ Several models have been proposed 
for the lattice dynamics of carbon nanotubes, ranging 
from zone folding and force-constant models f^^i^i^i^i^ 
valence force-field models^i^i^ and tight binding^ to 



ah initio calculations fi^i^>2^i^ 

In this work we concentrate on the well-established 
fourth-nearest-neighbor force-constant model by Jishi et 
al.^ It was developed and optimized for graphene 
and subsequently also adapted to carbon nanotubes by 
Saito et al.^^ It has been reparametrized several times 
for graphene^^i^i^ but, to our knowledge, no further cal- 
culations for carbon nanotubes have been done. While 
the original parameters were empirically determined by 
fitting experimental data of graphite, we propose a 
new parametrization to fit the ah initio phonon disper- 
sion of graphene.— With some corrections, we use this 
parametrization also for the calculation of the phonon 
dispersion of achiral carbon nanotubes. Main atten- 
tion is paid to the long-wavelength acoustic modes and 
to the controversial question of the dispersion law of 
the transverse acoustic (TA) or flexure mode— While 
Saito et al^ obtain four linear-dispersing acoustic modes 
{uj cx g), we obtain two linear modes and a doubly de- 
generate quadratic dispersing flexure mode {oj cx q^). 
The quadratic dependence of the flexure modes of car- 
bon nanotubes, predicted by continuum models^Si^i^ 
and obtained by several ab initio calculations ^^i^i^ 
has been reproduced only by few force-constant mod- 
g^g^ 4o,43,52 Qj^jj q£ ^YiQ latter (the work of Mahan and 
Jeon^) pursued a detailed study of the symmetry rules 
that lead to a quadratic flexure mode and achieved the 
correct behavior with a three-parameter spring and mass 
model. It is chosen for comparisons throughout this pa- 
per. 

We concentrate on an accurate description of the 
acoustic phonons since they allow to predict and inter- 
pret several low-temperature thermodynamic properties 
and to prove quantum size effects in carbon nanotubes. 
Using our force-constant model we calculate the specific 
heat and the thermal conductance of carbon nanotubes of 
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different diameter and chirality. Although the exact dis- 
persion law of the acoustic modes is apparently irrelevant 
for the quantized thermal conductance,— the quadratic 
dispersion of the flexure modes results in a very different 
behavior of the low-temperature specific heati^ Several 
experimental and theoretical studies have been achieved 
for determining the specific heat of carbon nanotubes 
and nanotubc ropes. Most of the m^^i^^i^^i^"^ cover a wide 
temperature interval and provide approximate estima- 
tions of the power law of the T-dependence, but only 
few work o^^'^^ extend the temperature range down to 
~0.1 K and provide a more precise analysis. Popo\^ 
showed within a force-constant model that a T"'^ depen- 
dence of the specific heat at very low temperatures can 
be directly related to the flexure modes. Lasjaunias et 
al^ determined experimentally on a sample of nanotubc 
ropes that under ^^1 K dominates a dependence. 
We study the exact power law including temperatures in 
the mK region and illustrate how it is correlated with the 
acoustic modes and the dimensionality of the system. As 
in Ref. we prove the ID quantization of the phonon 
subbands in nanotubes. 

Quantum size effects in carbon nanotubes are ob- 
served even in the thermal conductance. The thermal 
conductance of phonon waveguides in the ballistic, one- 
dimensional limit has been calculated by Rego et al^ 
using the Landaucr formula and has been proved experi- 
mentally by Schwab et al..— Within the same formalism, 
we show that the phonon thermal conductance of carbon 
nanotubes is quantized and determine the thermal con- 
ductance quantum. Our results are in good agreement 
with theoretical calculations . ^"^'^^ 

This paper is structured as follows. In Sec. II we pro- 
vide a brief description of the system. Section III deals 
entirely with vibrational properties: After a summary of 
lattice dynamics and the model approach in Sec. Ill A- 
B, we present our results for the phonon dispersions of 
graphene and carbon nanotubes in Sec. Ill C-D. In Sec. 
IV we consider thermal properties: The basic concepts 
are recalled in Sec. IV A, while Sec. IV B-C show our 
calculations of the low-temperature specific heat and the 
thermal conductance, respectively. Section V contains 
final remarks. 



II. SYSTEM 

The typical honeycomb structure of graphene is de- 
fined by a 2D hexagonal lattice with a basis of two atoms, 
which we call atom A and B. The lattice vectors are given 
by ai = (\/3a/2, a/2) and as = (V3a/2, -a/2) with the 
lattice constant a = 2.46 A, as illustrated in Fig. [T] 

A carbon nanotube can be thought of as a single 
graphene sheet that is wrapped into a seamless cylin- 
der. It is common to define a circumferential vector 
and a vector parallel to the tube axisi^i^ The first one, 
called chiral vector, is defined in terms of the unit vec- 
tors of graphene, C — na.i + ma.2, and the sheet is rolled 



up in such a way that it becomes the circumference of 
the tube. The pair of integers (n, m) uniquely defines 
a particular nanotube and thus provides a classification 
among nanotubes. The translational vector T is per- 
pendicular to C and reproduces the periodicity of the 
nanotube structure along the axis direction. It reads 
T = ((2m 4- n)/dji) ai — ((2n -|- m)/dii) as, being da the 
greatest common divisor of {2m + n) and (2n + m), as 
T should be the smallest lattice vector in its direction. 
Alternatively a nanotube can be defined also by its ra- 
dius R and the chiral angle Q, which is given by the chi- 
ral vector measured relative to the direction defined by 
ai. In this work we concentrate on the particular cases 
of zigzag (6 — 0) and armchair (6 = n/6) nanotubes, 
which represent the class of achiral nanotubes. 

III. VIBRATIONAL PROPERTIES 

A. Lattice dynamics 

To derive the equations of motion for the atoms we 
use Hamiltonian mechanics, treating the atoms as point 
masses moving according to the laws of classical mechan- 
ics. To describe the ion configuration, characterized by 
the instantaneous location of the atoms, we use the fol- 
lowing notation for a crystal with a monoatomic basis 

R„(t) = R„ + u„(t). (1) 

Hence at time t the ion is located at R„(t), while R„ is its 
equilibrium position. In the limit of small displacements 
u„ of the atoms from their equilibrium position, the so- 
called harmonic approximation, the equations of motion 
are a set of coupled second order differential equations 
given by 

A/u„ = - ^ $(R„, R„) • u,„ (2) 

m 




armchair 



FIG. 1: The grapliene lioneycomb lattice with lattice vectors 
ai and a.2. A carbon nanotube can be constructed by rolling 
up the graphene sheet along C so that points O and A coincide 
(as well as B and B' do). O denotes the chiral angle. (Figure 
taken from Ref. [60l]l 



where M is the mass of the constituent atom and 
<i?(R„, Rm) is the 3x3 force- constant tensor that couples 
atom n and m. Due to lattice periodicity it is possible 
to search for solutions of Bloch-wave type 



(3) 



where A gives the amplitude of the mode, lo the fre- 
quency, and q the wave vector. Inserting Eq. ([3]) in ([2|). 
the equations of motion become 



w^A/A = ^$(R„,R,„)e"i-(^'"-^")A. 

m 

These can be written in the compact form 
D(q) • A ^ uj^A 



(4) 



(5) 



where we introduced the discrete Fourier transform 

m 

The matrix D(q) is called dynamical matrix and is a 
hermitian and positive definite matrix. In order to ob- 
tain the eigenvalues w^(q) and the eigenvectors A(q) we 
have to solve the secular problem detD(q) = for each 
q vector chosen according to the Born-von-Karman pe- 
riodic boundary conditions. The generalization of the 
above described theory to non-monoatomic basis systems 
is easy and can be found, e.g. , in Ref. [Gj]- Considering 
a three-dimensional system with r atoms per unit cell, 
the dynamical matrix has 3r x 3r components so that for 
each wave vector q there are 3r frequencies ujs{(\}, with 
s = 1, . . . , 3r. 



B. Force-constant model 

A practical method of investigating vibrational prop- 
erties of graphene and carbon nanotubes (CNTs) is given 
by phenomenological lattice-dynamical models. These 
try to construct the force-constant tensor starting by an 
analytic expression for the interaction energy of two or 
more carbon atomsj^SiSiiS or alternatively by approxi- 
mating directly the interatomic force constants by fitting 
experimental datai^i^ Such empirical models are based 
on a few adjustable parameters and are able to provide 
reliable information that is complementary to that ob- 
tainable from more advanced methods. Indeed, an alter- 
native tool is given by first principles or ah initio calcu- 
lations based on the quantum mechanical description of 
electrons i^^i^ This method does not rely on input from 
experimental informations and includes all relevant ef- 
fects, providing accurate, experimentally-confirmed and 
therefore ve ry p redictive results, as has been shown, e.g. , 
However, the computational effort is 



Ref. 



large, leading to several restrictions in particular for com- 
plex systems of considerable size. The advantage of the 
phenomenological models consists in their simplicity and 




FIG. 2: An atom A and its first nearest-neighbor atoms Bp 
{p = 1,2,3). (pr, 4>ti, and ipto represent forces in radial, in- 
plane, and out-of-plane direction. 



the possibility of fast application to almost every system. 
In view of the aim of our work, force-constant models 
turn out as the best choice for two reasons: (i) they pro- 
vide quick and reliable implementation for several CNTs 
of different diameter and chirality, (ii) they reproduce 
with a high level of accuracy especially the low-energy 
acoustic modes, which in turn determine almost entirely 
the low-temperature thermal properties of CNTs. 

We calculate the phonon modes of graphene and 
carbon nanotubes using the force-constant model pro- 
posed by Saito et al.^ This model consists in the di- 
rect parametrization of the diagonal real-space force con- 
stants including up to fourth nearest-neighbor interac- 
tions (4NNFC approach). This leads to a set of twelve 
adjustable parameters. The truncation after the fourth 
nearest neighbors is justified by the rapid decay of the 
force constants i"^^'^^ In the 4NNFC approach the force- 
constant tensor describing the interaction between an 
atom and its nth nearest neighbor on an arbitrarily cho- 
sen axis {e.g. , the x axis) has diagonal form 



<P = 




(7) 



where and (/)("'' represent the force-constant 

parameters in the radial (bond-stretching), in-plane, and 
out-of-plane tangential (bond-bending) directions of the 
nth nearest neighbors. The radial direction corresponds 
to the direction of the bonds and the two tangential di- 
rections are perpendicular to it, as illustrated in Fig. [21 
The force-constant tensors for nearest-neighbor atoms of 
the same neighbor shell, that are not located on the x 
axis, can be obtained by unitary rotation of the ten- 
sor of Eq. (O. The formalism is described accurately 
in Ref. [38|. 

For example for first nearest neighbors (n — 1) we 
obtain the force-constant tensor ^^^'^p^ between atom 
A and its neighbor Bp (p = 2, 3) by 

^(A,Bp) ^ u^^iOp) $(^-^1) U, [Op) (8) 

where Uz{6p) is a unitary rotation matrix around the z 
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sin(6lp) 
cos{9p) 




1 



(9) 



and 6p is the angle defined by atom Bl, A, and Bp. 

The above mentioned tensors describe the interaction 
between atoms in the plane and for carbon nanotubes 
these have to be adapted because of the curvature of 
the walls. It is possible to generate the force-constant 
tensors for the atoms of the nanotube unit cell by rota- 
tion of the chemical bond from the two-dimensional plane 
of graphene to the three-dimensional coordinates of the 
nanotubci^ It is possible to generate the force-constant 
tensors for all the atoms of the unit cell of a nanotube 
from those related to one single atom of the cell, e.g. atom 
Al of Fig. [21 For atoms of type A the tensors of atom 
Al must be rotated by an angle '^i around the axis of 
the nanotube (here the y axis) 



(10) 



is the polar angle between Al and Ai around the 
circumference, li (p — i + I) is negative or zero we use 

+ p — i + I) instead of it. For atoms of type B we 
must rotate the tensors of atom Al first by tt around the 
z axis, and then by ^'^ around the y axis, as before 

<^(B',Ap) ^ U~\Tr) ^iAl,Bp~,+l) ^^^^^ 

(11) 

where Uy{i})) and Uz{ip) are unitary rotation matrices 
around the y- and z axis, analogously to that of Eq. 
The dynamical matrix is obtained by multiplying the 
force constant tensors obtained above by exjp{iqznT), 
where n is the number of the unit cell in which atom 
Al is situated, and T ~ |T| is the modulus of the trans- 
lational vector. 



C. Results for graphene 

In the calculation of the phonon-dispersion relation 
of graphene done by Saito et al.,— the force-constant 
parameters were empirically determined by fitting ex- 
perimental data of graphite obtained by inelastic neu- 
tron scattering. We perform, instead, a parameter fit to 
the ab initio dispersion relation of graphene calculated 
within density-functional perturbation theory by Bohnen 
and IIeid42, The corresponding sets of force constants are 
listed in Table U 

Precisely, the fitting process was performed through a 
matching of the force-constant tensor to those obtained 
by ab initio calculations 4^ This procedure is however 
limited by the constraint of including up to fourth near- 
est neighbors in the force-constant tensor, as required 
by the model approach. We varied and optimized the 
force constants in order to fit as closely as possible the ab 
initio phonon dispersion. We follow Gartstein^ choos- 
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FIG. 3: Dispersion relation of graphene calculated with the 
fourth-nearest-neighbor model (4NNFC) (solid lines) in di- 
rect comparison with ab initio calculations (dotted lines) of 
Bohnen and Held.— (a) 4NNFC approach with the original 
parametrization of Saito et al; (b) 4NNFC approach with 
our parametrization. The corresponding sets of parameters 
are listed in Table HI 



ing the in- and out-of-plane tangential force constants 



so as to satisfy 



6p) 



(4) 



0. 



This equality is required by the rotational invariance of 
the graphene plane, and the original parameters of Saito 
and coworkers do not obey this rule. Figure [3] shows the 
phonon-dispersion relation resulting from: (a) the orig- 
inal parametrization and (b) our new parametrization. 
Both curves are superposed to the ab initio dispersion of 
Bohnen and Held (dotted lines), for direct comparison. 

The dispersion relation of graphene comprises three 
acoustic (A) and three optical (O) modes, which are ei- 
ther out-of-plane (Z), in-plane longitudinal (L) or trans- 
verse (T). The acoustic ZA mode shows a energy dis- 
persion near F rather than the linear dispersion of the 
TA and LA mode, which is typical for acoustic modes 
in 3D solids. The quadratic dispersion is a characteris- 
tic feature of the phonon dispersion of layered crystals^ 
and can be explained as a consequence of the Dgft point- 
group symmetry of graphenci^ Another consequence of 
the symmetry are the linear crossings of the ZA/ZO 
modes and the LA/LO modes at the K point. With 
respect to the phonon dispersion obtained by Saito et 
al. and in comparison with first-principles results, our 
new parametrization yields a considerable improvement 
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Neighbor shell 


Parameters by Saito et al. (Ref. [3211 

-^1"' 0^:' 


(pr 


Our parametrization 

<Pti Vto 


1st 


36.50 


24.50 


9.82 


41.8 


15.2 10.2 


2nd 


8.80 


-3.23 


-0.40 


7.6 


-4.35 -1.08 


3rd 


3.00 


-5.25 


0.15 


-0.15 


3.39 1.0 


4th 


-1.92 


2.29 


-0.58 


-0.69 


-0.19 -0.55 



TABLE I: Force-constant parameters for graphene in units of 10'' dyn/cm — 10 N/m. 





Mode 


Ref. [49] 


This work 


Ref. [37] 






Ab initio LDA 




Ab initio GGA 
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ZO 


914 


925 


881 




LO/TO 


1576 


1583 


1554 


M 


ZA 


474 


469 


471 




TA 


632 


626 


626 




ZO 


641 


633 


635 




LA 


1351 


1315 


1328 




LO 


1375 


1351 


1340 




TO 


1437 


1341 


1390 


K 


ZA/ZO 


538 


539 


535 




TA 


1010 


969 


997 




LA/LO 


1243 


1208 


1213 




TO 


1302 


1408 


1288 



TABLE II: Phonon frequencies (cm ') of graphene at high 
symmetry points. 



in the overall phonon dispersion. In particular the acous- 
tic modes provide a remarkable good fit to ab initio data. 
The frequency values at high symmetry points F, M and 
K (listed in Table [n| differ only by up to 4% from ab 
initio data, with exception of the TO mode (6.7% at 
M and 8% at K), which will be discussed later. In the 
high-energy range, our parametrization leads to a quali- 
tatively correct rearrangement of the LA and LO modes 
along the line M-K and an improvement concerning the 
crossing of the LO and TO branches along the F-M and 
F-K directions. Nevertheless, there are still major di- 
vergencies from the ab initio dispersion for the LO and 
TO mode. Neither parametrization reproduces the initial 
upward curvature, called overbending, of the LO branch 
away from F that is observed in both the experimental 
dat a^^'^^i^°'^^ and in all published first-principles calcu- 
lations Furthermore, the TO phonon at 
the K-point is significantly higher than in the ab initio 
dispersion. It is known from literature that the highest 
optical phonon branch is shaped by the effect of electron- 
phonon interactions, which results in a discontinuity in 
the frequency derivative at F and Kfi^ These discontinu- 
ities are called Kohn anomalies and are revealed by two 
sharp kinks in the phonon dispersion. The two Kohn 
anomalies originate from a non-analytic behavior of the 
phonon dispersion, which is impossible to be reproduced 
by a finite set of interatomic force constants. All few- 
nearest-neighbor force-constant approaches yield a con- 
tinuous slope at F and K. 

In summary, the divergencies from ab initio curves ap- 




0.5 LO 0.5 I.O 

qTlTl qTlTl 



FIG. 4: Phonon dispersion for a (10,10) CNT calculated with 
the fourth-nearest-neighbor model with: (a) The parametriza- 
tion of Saito et al.; (b) Our new parametrization. The corre- 
sponding parameters are listed in Table H] Both parametriza- 
tions have been subsequently corrected in order to obtain 
ljj = Q &t q = Q for the acoustic TW mode. 



pearing in the high-frequency region are due to the natu- 
ral limit of accuracy of empirical force-constant models, 
which consider a finite number of nearest-neighbor atoms 
and miss the long-range character of the dynamical ma- 
trix. 



D. Results for carbon nanotubes 

In this section we present phonon dispersion rela- 
tions of achiral carbon nanotubes that rely on the ear- 
lier determined force constants of graphene. We con- 
centrate in particular on the commonly studied (10,10) 
and (10,0) CNT. The calculated phonon spectrum of 
a (10,10) CNT is illustrated in Fig. H for both the 
original parametrization of Saito et al^ and our new 
parametrization. In Ref. (ssj was proposed a scaling for 
the force-constant parameters in order to treat the curva- 
ture effect when rolling up the graphene sheet to form a 
nanotube. This scaling prevents from obtaining a wrong 
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shift of the rotational acoustic mode (TW) at g = from 
about « 4 cm~^. In the present work we do not apply 
the same rescaling, but vary empirically only the out-of- 
plane tangential force constants ^l""* = 4>)!^ (1 + e*^"^) of 
graphene. These are responsible for vibrations perpen- 
dicular to the atom-bonding plane and thus are the most 
subjected to changes when rolling a plane sheet into a 
cylinder. This effect increases with decreasing tube diam- 
eter. For a (10,10) nanotube and Saito's parametrization 
we obtain a frequency of the twisting mode of w ~ 10"'^ 
cm~^ at (J = by varying only cfif^ ili^ At this point it 
is important to observe that even very small variations 
in the force constants can have considerable effects on 
the low-frequency modes. In particular the frequency of 
the quadratically-dispersing modes are strongly affected 
by modifications of the parameters and can even become 
imaginary. For this reason, with our parameter set it 
was not sufficient to correct (fito only for the fourth near- 
est neighbors, but also for third and second neighborsi^ 
For the latter the correction is smaller, because these are 
less affected by the effect of curvature. It results in a 
frequency of the TW mode lo ~ 10~^ cm~^ at q = 0. 

The condition of infinitesimal rotational- and transla- 
tional invariance imposed on the force-constant tensors 
gives rise to four zero-frequency modes at q = 0. Near 
the F point the highest-energy acoustic mode is the lon- 
gitudinal (LA) mode, followed, in order, by a twisting or 
torsional mode (TW) and a doubly degenerate transverse 
or flexure mode (TA) . Figure [5] shows in detail the low- 
energy region of the phonon spectrum of a (10,10) CNT 
for three different cases: (a) and (b) are calculated with 
the 4NNFC model with the parametrization for graphene 
of Saito et al^ and with our parametrization, respec- 
tively, both corrected for nanotubes. Panel (c) is calcu- 
lated with a three-parameter spring-and-mass model for 
carbon nanotubes presented by Mahan and Jeon^ and 
is shown for direct comparison. The 4NNFC model with 
the original parametrization of Saito et al. shows a lin- 
ear dispersion at small wave vectors for all four acoustic 
modes. With our new parametrization of the 4NNFC 
model, adapted to nanotubes, we found that, while the 
high-energy optical phonons do not vary significantly, re- 
markable changes occur for the acoustic modes. The 
two degenerate TA modes now show quadratic disper- 
sion near the zone center, which was not given by Saito's 
parametrization. The model of Mahan and Jeon also ob- 
tains the quadratic behavior of the flexure modes, due 
to a detailed analysis and implementation of symmetry 
rules, which arc required by the correct force constants. 

Also in case of nanotubes with other diameters and 
chiralities we obtain the correct quadratic dispersion 
when applying our new parametrization. Figure [S] shows 
the phonon dispersion of a (10,0) CNT calculated with 
the 4NNFC model: Saito's parametrization provides lin- 
ear dispersions for all four acoustic modes (panel (a)-(c)), 
while with our parametrization we obtain the quadratic 
TA mode (panel (b)-(d)), which is given also by the 
model of Mahan and Jcon (panel (e)). In particular. 
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FIG. 5: Low- frequency region of the phonon dispersion re- 
lation of a (10,10) CNT, shown near the F point, calculated 
within: (a) 4NNFC model with the original parametrization, 
(b) 4NNFC model with our new parametrization, (c) Force- 
constant model by Mahan and Jeon^ based on three free pa- 
rameters. While in (a) all four acoustic branches have linear 
dispersion for small wave vectors, in (b) and (c) two increase 
linearly with g (LA, TW), and two are degenerate (TA) and 
increase quadratically (TA). 



after adapting the force constants (/jI"-*^ as in the case 
of the (10,10) CNT, we obtain uj ~ lO'^ cm"! at g = 
for the TW acoustic mode, for both parametrizations. 

Furthermore, we concentrate on the important 
Raman-active radial breathing mode (RBM) . This mode 
arises from a radial expansion and contraction of the en- 
tire tube. It is unique to single-walled CNTs and plays 
an important role in experiments] ^^'^^ One of the most 
important applications of the RBM is the determina- 
tion of nanotubes diameters on the basis of Raman data, 
through the expected dependence of the RBM frequency 
on diameter 



(12) 



where Ci is a constant, C'2 possibly depends on the di- 
ameter dt and k is an exponent. This functional depen- 
dence was first introduced by Jishi et al^ with C2=0 
and K~l. Several articles and a range of values of Ci 
have been published, differing from each other by a few 
per cent. A review of the experimental and theoretical 
values can be found in Rcf. [61|. For isolated tubes the 
values ranges from Ci = 218 to 248 cm^^ nm. 

We verify the relation of Eq. and analyze the chi- 
rality dependence of the RBM frequency. For this pur- 
pose we calculate the RBM frequencies of a number of 
armchair and zigzag nanotubes with the 4NNFC model 
with our parametrization adapted to CNTs. The ob- 
tained frequencies are almost perfectly inverse propor- 
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FIG. 6: Phonon dispersion for a (10,0) CNT: (a) and (c) 
4NNFC model with the parametrization of Saito et al; (b) 
and (d) 4NNFC model with our new parametrization. The 
two parametrizations show respectively linear and quadratic 
dispersion for the TA mode. For both the 4>to constants have 
been corrected in order to obtain zero frequency for the acous- 
tic TW mode at the Brillouin-zone center, (e) Low-frequency 
region of the spectrum, calculated with the model of Mahan 
and Jeon, which gets the correct quadratic dispersion of the 
TA mode. 



tional to the radius of the tube, as shown in Fig. [7^, and 
independent on chirality. The RBM frequency decreases 
with increasing tube diameter and becomes zero in the 
hmit of infinite diameter, which corresponds to the out- 
of-plane tangential acoustic mode of graphene at g = 0. 
By fitting the frequencies of the RBM to tube diameters 
by the relation of Eq. (|12|1. we get neghgible values for 
C2 (order of 10~^ cm~^ nm), and Ci = 212 cm^^ nm for 
armchair tubes (Fig. [7j3), which is in satisfactory agree- 
ment with the experimental value of 224 cm~^ nm.— For 
zigzag nanotubes (Fig. [71;) we obtain Ci = 209 cm~^ nm. 
The values of Ci are in agreement also with previous cal- 



FIG. 7: (a) Frequency of the radial-breathing mode of various 
armchair (n = 3 — 12) and zigzag (n = 6 — 20) tubes as a func- 
tion of the nanotube diameter, calculated with the 4NNFC 
model with our new parametrization. For comparison, the 
dashed lines show the results for the original parametrization 
of Saito. (b) and (c) Frequency of the RBM as a function 
of the inverse tube diameter for armchair and zigzag tubes, 
respectively. The solid lines are a linear fit to the data ex- 
cluding the small-diameter tubes (3,3), (4,4), (6,0), and (7,0), 
which are marked by circles. These show a deviation from the 
predicted behavior, with a decrease in the RBM frequency. 
Reference [sj explains it as a consequence of the hybridiza- 
tion changes and the decrease of the tt interaction induced by 
the curvature. 



culations . ' 1 The frequencies obtained by the original 
parametrization of Saito et al. are displayed in Fig.[7|by 
the dashed lines. We obtain Ci = 223 cm~^ nm for arm- 
chair tubes and Ci = 222 cm~^ nm for zigzag nanotubes. 
They arc in better agreement with experimental values 
than our parametrization, since the RBM frequencies are 
on average about 8cm^^ higher. 

A possible chirality dependence is below the resolution 
of the data. Indeed the proportionality constant Ci dif- 
fers only by about 1% between armchair and zigzag nan- 
otubes. This can be explained by the fact that the RBM 
corresponds to a stretching of the graphene sheet in the 
[110] (armchair tubes) or [100] (zigzag tubes) direction. 
Because the system is isotropic in the hexagonal plane, 
the elastic constant that the describes the stretching of 
a graphene sheet is independent on the direction.— 

Beneath the RBM, also the other low-frequency modes 
depend strongly on the tube diameter.— Instead, the 
higher-frequency modes do not have such a strong di- 
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ameter dependence since their frequencies are more sen- 
sitively determined by the local displacements of the 
atoms. 



IV. LOW-TEMPERATURE THERMAL 
PROPERTIES 

A. Methods 

1. Specific heat 

In order to characterize the specific heat at constant 
volume of low-diameter single wall CNTs, we start from 
the definition 



cv 



1 f dE 



V \dT 



(13) 



where the internal energy E ~ k-^T^ {dhiZ / dT) is defined 
through the vibrational partition function Z . According 
to statistical thermodynamics, the partition function of 
a system of independent harmonic oscillators can be di- 
rectly expressed in terms of the phonon frequencies by 



nr 



-hL0s{q)/2kTiT 
Q-huj,{ci)/kBT ■ 



(14) 



Here, q is the phonon wave vector. Ws(q) are the phonon 
frequencies with mode index s = 1, . . . , 3r (r is the num- 
ber of atoms per unit cell, see Sec. IIII A[) . T is the tem- 
perature, /cb the Boltzmann constant and h the Planck 
constant. The specific heat can thus be written as 



Cv 



E 



fti^s (q) 

2kBT J smh'^{huJs{q)/2kBT)' 



(15) 



We do not distinguish between specific heat at constant 
volume, Cv, or constant pressure, cp, since the model ap- 
proaches do not include thermal expansion of the lattice. 
Anyway, the difference between cy and cp is only in the 
range of a few percent]^ For densely spaced values of 
the wave vector q it is possible to replace the sum by an 
integral 



V 



(27r)3 
3rN I g{u})du} 



(16) 



(17) 



where g{uj) is the density of states (DOS) and N is the 
number of unit cells. The expression for the specific heat 
results 



Cy = Srfcj 



duj 



2kBTj sinh\nw/2kBT) 



(18) 



Therefore the specific heat depends in a detailed way on 
the frequency spectrum g{ijj) of the normal modes. 



2. Landauer phonon transport 

Phonon heat transport in mesoscopic systems can be 
investigated using methods analogous to the Landauer 
description of electrical conductancc i^^'^^ We consider a 
model of an ideal one-dimensional heat conductor, built 
by two long perfect leads that join a central segment in 
which the phonon scattering occurs. Only elastic scat- 
tering is taken into account, while phonon-phonon inter- 
action is neglected. The free ends of the two leads are 
connected to reservoirs of temperature Thot and Tcoid, 
respectively. No scattering occurs at the reservoir-lead 
connections. The energy flux of the right/left moving 
phonons is given by^^i^liZ^ 

1 f°° 

„ Jo cold 



2tt 



^ / dujhi^riY,at/iuj)%{uj) (19) 



cold 



where LUs{q) is the dispersion relation of the discrete 
mode s, Vs{q) is the group velocity and Ts{q) are trans- 
mission coefficients characterizing the coupling of wave- 
guide modes to the reservoirs. The total heat current is 
therefore Jph = J"*" — J~ ■ Assuming perfectly adiabatic 
contact between the thermal reservoirs and the ballistic 
quantum wire, the transmission function for a monoton- 
ically dispersing mode s is the step function 



for L 
otherwise. 



< w < 



(20) 



Instead, for non-monotonic dispersions, given a frequency 
U, the transmission 7l(Zu) is defined as the number of 
crossings of the line uj ~ lJ with the phonon dispersion of 
the mode s. With the total transmission function given 
by T = 7^(a;), the Landauer energy flux resultsS^ 



J, 



ph 



duj 
2^ 



Eventually, the thermal conductance is defined as 



Jph 

AT 



(21) 



(22) 



with AT = Thot — Tcoid- In the limit of linear response, 
AT < T = (Thot + Tcoid)/2, we obtain using Eq. ^ 
and the substitution x = fiw/k-BT 



'^ph 



fc|T 
h 



dx ■ 



x'^c^ 
(e^ - 1)2 



■ rix 



knT 



(23) 



This equation plays the role of a 'universal' phonon con- 
ductance. An important statement is that the result is 
independent of all details of the dispersion curve except 
the transmission function. This arises because the den- 
sity of states in the frequency integral is canceled by the 
group velocity. 
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B. Specific heat results 

The specific heat of carbon nanotubcs is mainly de- 
termined by phonons, while electronic contributions to 
it can be neglected even at a few Kelvini^i According to 
Eq. it depends sensitively on the characteristics of 
the phonon spectrum and on its vibrational density of 
states (DOS). 

The specific heat calculated from the theoretical DOS 
spectra is shown in Fig. [5] as a function of temperature. 
At high temperatures, the specific heat of all the differ- 
ent approaches and all chiralities converges to the classi- 
cal limit of Sfce/Af = 2078 mJ/gK with M being the 
atomic mass of carbon (see inset in Fig. [S^). In the 
low-temperature regime that we are interested in (be- 
low 20 K. sec Fig. [8)d); the specific heat of graphenc is 
dominated by the quadratic out-of-plane bending mode 
and is expected to have a linear T dependence at very 
low temperatures. The nanotube curve is lower than 
the graphene one because the tube has no low-energy 
counterpart to the layer-bending modesi^ In this tem- 
perature range, nanotube modes with ^^^^(q) ^ UbT 
will negligibly contribute to Eq. (fT8|) , since the integrand 
will vanish exponentially. Optical phonons are frozen out 
and only the long-wavelength acoustic modes are pop- 
ulated, because for these holds ijJs{q) ^ as g ^ 0. 
Therefore, the acoustic modes alone determine the low- 
temperature behavior of the specific heat. Only at a tem- 
perature Topt ~ ^kJopt/Gfee does the lowest- lying optical 
subband with frequency Wopt at g = begin to contribute 
to the specific hcatj^ The two cy curves calculated by 
the model of Mahan (dot-dashed line in Fig. [SJd) and 
by the 4NNFC model with our parametrization (solid 
line) coincide in the temperature range of a few Kelvin. 
This was expected because both models predict quadrat- 
ically dispersing flexure modes and thus analogous low- 
temperature behavior. The slope of the curves increases 
smoothly when the first optical subband begins to con- 
tribute to Cy. In the model of Mahan this takes place 
approximately at Topt ~3 K, while for our parametriza- 
tion it appears about 5 K, due to the different frequencies 
of Wopt, that are 12.6 cm~^ and 20.2 cm~^, respectively. 
The 4NNFC model with Saito's parametrization instead 
yields linearly dispersing flexure modes and the (logcy) 
vs. (logT) curve (green line in Fig. (Bja) shows a higher 
slope than the two curves described just now. Due to 
the lowest lying optical mode with Wopt = 21.0 cm^^, 
the slope increases at ^5 K, as expected. 

The low-temperature behavior of the cy vs. T curves 
needs to be analyzed in more detail. Indeed, the behavior 
of Cy contains informations regarding the dimensionality 
of the system through a fixed correlation between cy, the 
DOS and the exact dispersion lawi^ According to Ta- 
blc lllll through evaluation of the exponent a in the power 
law Cy a it is possible to get informations about the 
dimensionality of the system. Since nanotubes are quasi- 
one dimensional (ID) systems consisting of roUcd-up 2D 
sheets, they should display both ID quantum size effects 
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FIG. 8; (a) Specific heat as a function of temperature for 
graphene (dotted line), calculated with the 4NNFC model 
and our new parameters, and for a (10,10) CNT calculated 
with the model of Mahan- Jeon (dot-dashed Une), the 4NNFC 
model with our parameters (solid line), and the 4NNFC model 
with the original constants of Saito et al. (dashed line). The 
inset shows a wider temperature interval; cy approaches the 
value 2078 mJ/gK for high temperatures, (b) Specific heat 
on a logarithmic scale for the low frequency region. The as- 
signment of the line styles is the same as in (a) . 



Dimension 


Phonon 


Phonon 


Specific 




dispersion 


DOS 


heat 


ID 


LJ cc 


g{uj) cc 


cv cc ^/^" 




uj cc q 


g{uj) = const 


cv cc T 


2D 


U! cc q'' 


g{i^) = const 


Cv cc T 




u) cc q 


g{uj) cc ui 


Cv oc 


3D 


U! (X q 


g{ijj) oc oj'' 


Cv oc 



TABLE III: Low-temperature behavior of the specific heat. 
The dimensionality of the system is correlated to the density 
of states and, therefore, to the specific heat. At low temper- 
ature only acoustic modes are excited. These can have either 
linear or quadratic dispersion. 



and 2D features. 

Figure [5] shows the low-temperature specific heat for 
the (10,10) and the (10,0) with the 4NNFC model and 
our parametrization, and the inset shows the slope a 
of (logcy) vs. (logT) curves. In the mK temperature 
range a clearly tends to the value 1/2. This is almost 
entirely due to the doubly degenerate flexure mode with 
quadratic dispersion law (sec Table [ITT)) . which is domi- 
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FIG. 9: Low-temperature specific heat for a (10,10) and 
a (10,0) CNT calculated with the 4NNFC and our new 
parametrization. The inset shows the value of their slope 
a = d(logcv)/d(logr). 



nating in this temperature range. With increasing tem- 
perature (0.8 to 5 K for the (10,10) CNT), also the con- 
tribution of the two linearly dispersing modes becomes 
stronger and a holds values between 1/2 and 1 due to the 
superposition of these four modes. Above Topt the slope 
changes considerably due to the optical phonons and the 
tube is essentially 2D. This behavior is in accordance 
with theoretical predictions^ and is a direct confirmation 
of quantized ID phonon subbands in carbon nanotubes. 
Our results for low-temperature cy agree very well with 
previous calculations of Popov )^ and additionally we ex- 
tended the low-temperature limit by two orders of mag- 
nitude. We find good agreement also with experimental 
results of Lasjaunias et al.^ who measured the specific 
heat down to 0.1 K and fitted their measured curves with 
a power law of 0M3T°-^^ + O.OSST^. However, experi- 
mental measurements are usually performed on bundles 
of nanotubes, whose properties can differ greatly from 
those of isolated tubes. The adding of tubes to a bun- 
dle suppresses in particular the bending flexure modes, 
with a consequent increase of the exponent a in favor of 
a linear T dependence. 

We achieved similar results for other armchair and 
zigzag nanotubes. However, since the first optical sub- 
band edge varies from tube to tube and depends on the 
model, the turning points in cy vs. T curves are dif- 
ferent, causing a crossover of cy curves. The general 
uptrend and the high-temperature limit are the same. 
Figure [TUk shows the specific heat curves for a (10,10) 
and a (10,0) CNT. The tube diameter influences the spe- 
cific heat of carbon nanotubes, especially in the range of 
25-350 K. In order to determine the effect of tube diame- 
ter on the specific heat, additional results for T = 300 K 
are displayed in Fig. [TOb . using the 4NNFC model with 
our parametrization. At a fixed temperature the specific 
heat increases with increasing tube diameter. This was as 
expected, since for very large diameters the curve should 
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FIG. 10: Chirality dependence of the specific heat: (a) Tem- 
perature dependence of cv for a (10,10) CNT (straight lines) 
and a (10,0) CNT (dashed lines). The index 1 refers to the 
model of Mahan-Jeon and the index 2 to the 4NNFC model 
with our parameters; (b) The effect of tube diameter on zigzag 
and armchair CNT specific heat. At a given temperature the 
specific heat increases with the increase of tube diameter. The 
upper limit is given by graphene, with cv = 794 mJ/gK at 
300 K. 



approach the cy value of graphene, which is 794 mJ/gK 
at 300 K. However, the effect decreases at large tube 
diameter. The chirality shows only a small effect in the 
tubes specific heat, with cy of the zigzag tubes lying over 
that of the armchair tubes. This small effect is negligible 
and could also be caused by inaccuracies of the model de- 
scription. The results are in good agreement with those 

ofRcf. [HIES. 



C. Thermal conductance 

In the following, we demonstrate that at low temper- 
atures a carbon nanotube behaves as a ballistic, one- 
dimensional wire and that the phonon thermal con- 
ductance is quantized. The thermal conductance can 
be calculated by evaluating Eq. In this expres- 

sion, the integrand is given by the product of two func- 
tions: the transmission function and a weight function 
x^e^/(e^ — 1)^. The former is related to the phonon 
spectral properties of the nanotube and the latter takes 
into account the effects due to temperature. Figure [TT] 
shows the transmission function T = J^s'^'^i^) ^'-'^ ^ 
(10,10) CNT. Some branches count doubly because of 
their degeneracy. For high temperatures, these two func- 
tions arc non-zero within the same range, which means 
that all the transmission modes contribute to the thermal 
conductance. Whereas, in the limit of low temperature, 
the broadness of the weight function is extended only to 
the low-energy modes. Therefore, at low temperatures 
only four acoustic modes give an appreciable contribu- 
tion to the thermal conductance of a carbon nanotube. 
In this temperature regime Eq. (j23p becomes greatly sim- 
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FIG. 11: Phonon dispersion of a (10,10) CNT and the trans- 
mission function T = '^^Ts{ui). The latter is a sum of s = 120 
step functions. 



plificd 



klT 



dx 



(c^ - 1)^ 



= 4. 



3h 



(24) 



Here, the factor 4 represents the number of acoustic 
modes. The upper limit of the integral is of few impor- 
tance, because the integrand function falls off rapidly, 
before the successive step in the transmission function 
takes place. From Eq. (pi)) results that a fundamental 
relation holds for each mode 



3h 



(25) 



This quantum of thermal conductance represents the 
maximum possible value of energy transported per phonon 
mode. It docs not depend on particle statistics, therefore, 
is universal for fcrmions, bosons, and anyons/i^ Further- 
more, it is independent of any material parameters and 
of precise details of the dispersion law. This is clear since 
to construct the transmission function as in Eq. (PO)) . it 
does not matter whether the dispersions are linear or 
quadratic, but the branch upper and lower limits should 
be accurately computed. 

The phonon thermal conductance Kph of a 
(10,10) CNT as a function of temperature is shown in 
Fig. 112b ,. already normalized by Akq. For temperatures 
in the range of a few Kelvin, this ratio reaches the 
constant value of 1, independent of the model approach 
and of chirality. This behavior confirms that the thermal 
conductance of carbon nanotubes is quantized. Despite 
the quantization, the curves do not present steps because 
of the broadening of the Bose-Einstcin distribution in 
comparison with the energy gap between subband edges. 




y*^" (10,10) CNT 

— • Mahan model 

4NNFC model 

(parabolic TA mode) 

_ _ 4NNFC model 
(linear TA mode) 



5 10 15 20 
Temperature (K) 



25 



1.4 



1.3 



y° 1.2 
1.1 



1.0 



0.9 



' 1 ' 

: (b) 


(10,10)/ ush/^ y^- 


















1 . 1 . 1 . 1 . 1 . 



10 15 20 25 
Temperature (K) 



30 



FIG. 12: (a) Phonon thermal conductance for a (10,10) CNT 
calculated with the Mahan- Jeon model (dot-dashed line), 
the 4NNFC model with our parameters (solid line), and the 
4NNFC model with the original constants of Saito et al. 
(dashed line). The two latter show, respectively, quadratic 
and linear dispersion for the TA mode, (b) Thermal con- 
ductance for several carbon nanotubes, calculated with our 
parametrization. 



The length of the plateau depends on the lowest optical 
frequency of the dispersion curve. Indeed, the turning 
point in the curve of the (10,10) CNT calculated with 
the model of Mahan (Fig. \T% .. dot-dashed line) is about 
2K, while it is higher for the 4NNFC model with both 
paramctrizations (about 3K). As predicted, the exact 
dispersion law of the acoustic TA mode does not affect 
the qualitative behavior of the Kph vs. T curve. The 
two curves for the 4NNFC model result respectively 
from a quadratic and linear dispersion of the TA mode, 
but show the same low-temperature behavior. The 
deviations above 10 K are due to small differences in the 
respective optical frequencies. 

The results of the 4NNFC model arc believed to be 
more accurate than the ones obtained using Mahan and 
Jeon's method for the phonon dispersion, because the 
latter does not correctly reproduce the graphene disper- 
sion. Eventually, the thermal conductance depends only 
on the tube radius and not on chirality. Results for arm- 
chair tubes are very similar to the ones for zigzag tubes. 
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when same diameters are compared. This arises because 
the energy hw of the lowest-lying optical modes is de- 
termined only by the tube radius and decreases approx- 
imately according to ~ 1 /i?^ (see Ref. [s^ ) . Figure [T^ 
shows the thermal conductance for some armchair and 
zigzag nanotubes. 

Our results are in very good agreement with those 
of Yamamoto et al^ and Mingo et al.^ Experimental 
studies were achieved by Schwab et al.^ who observed 
the quantum thermal conductance in a nanofabricated 
ID structure, which behaves essentially like a phonon 
waveguide. 

Finally, we would like to express a word of caution 
to specify the limits of the Landauer model of heat con- 
duction. It must be stated that it describes an ideal- 
ized case of ballistic transport through a one-dimensional 
waveguide, where the phonon transmission occurs with- 
out scattering by defects or scattering at the reservoir- 
lead connection. The waveguide and the reservoirs are 
coupled adiabatically and anharmonicity and phonon- 
phonon interactions are neglected. These conditions are 
fulfilled only at low temperatures, where the phonon 
mean free path is limited only by the size of the sys- 
tem and anharmonic terms are small compared with the 
harmonic part of the Hamiltonian. However, as far as the 
system size exceeds the mean free path, which strongly 
depends on temperature, scattering of phonons due to 
anharmonic terms of the interatomic potential begins to 
decrease the conductivity and the transport ceases to be 
ballistic. Indeed, anharmonic terms give rise to phenom- 
ena as finite phonon lifetimes and interaction between 
phonons. These determine strongly the transport proper- 
ties at higher temperatures and are responsible for finite 
thermal conductivity. 

V. CONCLUSIONS 

We presented a combined theoretical investigation of 
both vibrational and thermal properties of graphene and 
carbon nanotubes within a force-constant model. First 
we fitted the phonon dispersion of graphene to that ob- 
tained with ab initio calculations by Bohnen and Heid^S, 



and found reasonable agreement for the overall disper- 
sion and good agreement for the acoustic modes. The 
frequency values at high symmetry points F, M, and K lie 
close to those obtained by various first-principles calcu- 
lations (about 4%, with the exception of only one mode). 
Then we presented results for the phonon spectra of achi- 
ral carbon nanotubes and focused on the low-frequency 
region. The dispersion of the doubly-degenerate flexure 
mode shows lo oi behavior at long wavelengths, as 
predicted by several theoretical works. Particular atten- 
tion has been paid to the radial-breathing mode, with 
a detailed analysis of the frequency- and chirality de- 
pendence on the tubes' diameter. On the basis of the 
so-obtained phonon spectra, we calculated the specific 
heat and the thermal conductance of carbon nanotubes. 
The quadratic dispersion of the flexure modes leads to 
a \/T dependence of the specific heat at very low tem- 
peratures. This is a direct confirmation of the one- 
dimensional behavior of carbon nanotubes at low tem- 
perature. Concerning heat transport, we showed that 
nanotubes can conduct heat by ballistic phonon propaga- 
tion. At low temperatures the thermal conductance for 
a single phonon channel approaches a maximum value 
of kq = TT^fcgT/S/i, which is the universal quantum of 
thermal conductance. We showed that for nanotubes of 
different diameter and chirality, the thermal conductance 
reaches the value 4ko for T — > 0, where the factor 4 is 
due to the four acoustic modes of a nanotube. All our 
results are in very good agreement with theoretical and 
experimental data available in literature. 
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